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Abstract 

Employing  the  system  of  Integro-dlf ferentlal 
equations  which  govern  a  plasma  in  the  guiding 
center  limit,  stability  criteria  are  found  for 
infinite,  homogeneous,  colllslonless ,  guiding  center 
plasma.   A  one-dimensional  fundamental  solution  is 
derived  for  the  linearized  system,  and  a  well- 
posedness  theorem  is  proved  for  the  corresponding 
Initial  value  problem. 
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1.   Introduction 

The  setting  for  this  Investigation  Is  a  colllslonless , 
charge -neutral  plasma,  consisting  of  single  species  of  Ions 
and  electrons,  permeated  by  a  strong,  slowly  varying  magnetic 
field.   It  Is  assumed  that  the  mass-to-charge  ratio,  m/e ,  Is 
small  compared  to  unity  for  both  classes,  and  that  all 

particle  orbits  may  consequently  be  described  by  the  lowest 

12  3  4 
order  guiding  center  theory.  ' 

Accordingly,  the  orbits  consist  of  a  rapid  circular 

gyration  with  small  Larmor  period  about  a  magnetic  field 

line  on  which  is  superimposed  both  a  drift  velocity  which  is 

Independent  of  the  particle,  and  a  kinetic  velocity  parallel 

to  the  field.   The  high  frequency  gyration  induces  a  local 

symmetry  in  the  pressure  tensor  as  a  result  of  the  circular 

motion  orthogonal  to  the  field  lines.   As  a  result,  we  obtain, 

in  general,  an  anisotropic  stress  tensor  which  depends  on 

two  scalar  pressures  and  the  orientation  of  the  unit  vector, 

p,  in  the  direction  of  the  magnetic  field.   Expressed  Invari- 

antly,   the  elements  of  the  tensor  are  given  by-^ 

(1)      p.  .  -  p  p.p  .  +  p, (6.  .  -  p.p  .) 

where  p.  are  the  components  of  p,  and  5.  .  is  the  Kronecker 
delta . 

The  model  v/e  employ  is  given  below,  and  was  derived 
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6  7 
by  H.  Grad.  '  ' 

(2)      0  =  (^+  ^-  t"+^PJ)f  +  ^  MP'°^  ^ T ^-) 


a  /I  ,2 


+  *^  '  ^  ^^  "^    ^'^^  +  v(u-k)  ^  f-]  , 


+  , 


O)  0  ^  11+  V  x(B  X  u)  , 


_^> 


(4)       0  =  pB  X  ^  +  B  X  dlv(P^+  P  )  -  n^  B  K  (  (V  ^  B)><  B) 

at;  ^.Q 


In  these  equations,  all  quantities  are  time-averaged  over  a 

+ 
Larmor  period.   The  distribution  functions,  f~,  are  taken  as 

functions  of  the  magnetic  moment  per  unit  mass,  |x,  and  the 

kinetic  speed,  v,  parallel  to  the  magnetic  field  line  B. 

They  are  normalized  as  mass  densities  in  \i-v    space  in  our 

+      + 

treatment  (/fdudv  =  p)-   The  quantities  1P~  and  p   are  the 

pressure  tensors  and  mass  densities  respectively.   (Super- 
scripts "*"  and  '  will  refer  to  ion  and  electron  quantities.) 
The  macroscopic  flow  velocity  orthogonal  to  B  is  u  and  k 
represents  the  curvature  of  B.   The  quantities  a"  are  defined 
by 


(5)       a 


*  -     (m/e)+ 


(m/e)^  +  (m/e)' 


Equations  (2)  are  the  collislonless  Vlasov  equations  for  the 


guiding  centers  of  the  ions  and  electrons.   Equations  (3) 

are  Maxwell's  equations  for  the  time-evolution  of  the 

magnetic  field  in  which,  consistent  with  the  lowest  order 

guiding  center  theory,  the  electric  field  has  been  replaced 

by  (B  X  u).    The  last  equations  of  the  system  are  the 

magnetohydrodynamic  equations  for  the  components  of  velocity 

orthogonal  to  B.   In  these,  the  current  has  been  replaced 

by  7^  (V  X  B). 
^o 
In  this  paper,  a  space-filling  plasma  with  homogeneous 

steady  state  is  assumed  and  our  system  is  linearized  about 
it.   The  initial  value  problem  for  the  linearized  system  is 
then  investigated. 

The  first  two  sections  exhibit  the  linearization,  and 
a  uniqueness  theorem  is  proven  for  the  purpose  of  justify- 
ing the  results  of  the  following  section.   In  that  section, 
which  contains  the  main  results  of  the  paper,  stability 
is  investigated  via  Fourier  and  Laplace  transforms.   The 
limitation  of  the  theory  to  slowly  varying  fields  is  mani- 
fested in  a  nondispersive  relation  which  determines  the 
exponentially  growing  modes  of  the  system,  with  the  conse- 
quence that  the  existence  of  such  modes  implies  the 

Q 

non-well-posedness    of  the  initial  value  problem.   Follow- 
ing the  discussion  of  stability  some  results  in  wave 
propagation  are  derived  for  special  lim.iting  values  of 
parameters  of  the  problem.   In  the  last  section,  a  funda- 


5  - 


mental  solution  for  the  one -dimensional  plasma  is  derived 
for  arbitrary  magnetic  field  orientation.   Some  asymptotic 
properties  of  this  solution  are  exhibited.   In  the  appendix 
we  treat  some  of  the  integrals  encountered  in  the  paper. 

A  word  must  now  be  said  concerning  the  stability  results 

which  will  be  derived.   In  our  investigation  we  employ  a 

q 
method  similar  to  that  of  0.  Penrose,^  and  show  that  the 

guiding  center  theory  predicts  stability  for  wide  classes 
of  steady-state  distributions.   In  particular,  it  will  be 
shown  that  the  two  necessary  conditions  for  stability  which 
are  obtained  from  the  guiding  center  theory  extend  those 
obtained  from  a  purely  macroscopic  model  by  R.  Liist,   and 
then  are  in  agreement  with  those  derived  in  a  special  case 
by  H.  Furth  who  considered  a  particular  field  perturbation 
and  linearized  the  Vlasov  equation  about  a  two-temperature 
Maxwellian  distribution.   It  is  also  shown  that  the  afore- 
mentioned conditions  are  not  sufficient  for  stability,  but 
when  coupled  with  a  third  inequality  which  is  shown  to  be 
satisfied  by  distribution  functions  that  are  simply-peaked, 
the  new  set  is  sufficient. 


Acknowledgement:   The  author  wishes  to  take  this  opportunity 

to  express  his  indebtedness  to  Prof.  Harold  Grad  for  suggesting 

this  investigation,  and  to  Prof.  Harold  Weitzner  for  his  many 
helpful  observations  and  suggestions. 
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2.   Linearization 

We  now  linearize  equations  (2) -(4)  about  a  homogeneous 
steady  state.   Quantities  associated  with  this  state  will 
be  distinguished  by  a  "o"  subscript  or  superscript.   Thus, 
B   Is  the  magnetic  field  vector  In  the  steady-state  con- 
figuration.  The  absence  of  such  an  Index  on  a  dependent 
variable  will  Imply  a  time -dependent  quantity.   It  Is  to  be 
assumed  that  such  quantities  are  sufficiently  small 
(together  with  their  derivatives)  to  Insure  that  products 
of  them  may  be  neglected. 

Our  coordinate  system  Is  chosen  so  that  the  magnetic 

field  vector  B   Is  parallel  to  the  x-axls,  and  the  plasma 
o 

Is  at  rest  (i.e.,  u^  =  (0,0,0)).   Thus, 


B 


;b^  +  B^,  B^,  B3) 


(6) 


■u  =  (0,  u^,  u^) 

+    +    + 

P-  =  P^  +  9 

+    +        ±      ^ 

f-  =  f:(^i,v)  +  f  (^L,v;  x,t) 


To  first  order  we  have 


(7)      I  B  I  =  /(b^+B3_)2+b2  +  b|  =  B^  4-  B-^  . 
Therefore,  the  components  of  p"  =  b/1b|  are  given  to  first 


5  - 


order  by 


(8) 


P 


1 


^0  +  % 
Bo+% 


=  1  ; 


^2,3  _  ^2,3 


2,5  -  B^+B-L 


B. 


o 


Thus,  using  the  representation  of  the  pressure  tensor  as 
given  by  (l),  we  obtain  the  linearized  tensors 


(9)     P- 


+ 


,+ 


2 


+ 


p       +  /   V    f  ~di-Ldv 


I    o-         o^,    ^2 
(P^    -  P„     )    ^ 


+  ±      B^  ^+  B-i  + 

-(p°-  -  p°-)    ^        pr(l+  ^)+B Jnf-dHdv 


^■t    ^+   B-, 
/  o    o-N   5 

(Px  -P„  )  B~ 


0 


^ 


(Px  -P„  )b- 
o 


0 


r+  B. 

pr(i+B^) 

o 

+ 

+B  /  li-f  "dixdv 
o  •' 


In  deriving  (9)  we  have  used  the  following  integral  repre- 

12 

sentatlons  of  the  pressures 


>~ 


P 


dM. 


(10) 


II 

J 

0 

+ 

1  = 

= 

+ 


V  dv  f~  (m-,v) 


.+« 


B 


H-dM- 


J 
o 


+ 


dv  f   (m-,v) 


Substituting  these  expressions  in  our  system,  and  keeping 
only  first  order  terms,  we  obtain  the  linear  equations 
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(11) 


+ 


+ 


+ 


0  = 


+ 


+ 

^  +  V  4-)f"  +  ;  ([Tl]a~-  n)f"   -v^U-^ 
I  o  J 


^^  +  ^  ^^ 


+ 
o 


f:  ^Sr^   Sr, 


+  + 


—     +   ^W  a"  f"    ^ 
B   I  dXp   dx^y        o,v  dx-. 


2 
dM-v  dv 


\Po  p 


Oi 


(12) 


0 


0  = 


SB, 


\  dXr^        dx 


3 


5F~ 


^  2,3 


(13) 


Sr 


2'^   v^  .S  r. 


B 


+ 


Po^^ 


+  c 


-,+ 


2      ^^1 


3x 


2,3 


^L(f   +  f   )d^LdV 


in  which  we  have  used  the  abbreviations 


^2,3  -  ^o  "2,3 


(14) 


o+     0+ 


o-   o- 


(Px   -  Pn   )     (Px   -  Pn   ) 

— :+ + = — 


O 


B 


1   /  O  ,    Ox 

=  7^  (11-+  PJ 


^o    o 

b" 


1    /   O   ,     O        On 

-  (tt-  +  Px  -  Pj 


p 


^o   '^o 


o 


In  the  following  the  further  abbreviations 


(15) 


2/^+  .+ 

o 


<v^>=  //  dudv  v^(f+  p 
<(ia)  =  //  nd^Ldv(f"^  +  f") 


f"  P") 
^o^ 
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will  be  employed 


3-   Uniqueness.  . 

The  methods  of  the  following  section  rely  very  much 
on  a  Laplace  transform  analysis  of  equations  (11)-(13)- 
Such  an  analysis  is  limited  to  solutions  which  grow  no 
faster  than  some  exponential  which  is  linear  in  time. 

o 

However,  as  G.  Backus   has  noted  for  some  first  order 
differential  operators  in  time,  solutions  may  exist  which 
grow  faster  than  any  linear  exponential.   Thus,  a  unique- 
ness theorem  is  necessary  in  order  to  substantiate  the 
results  obtained  by  such  a  transform  analysis  for  the 
initial  value  problem. 

It  will  be  assumed  that  initial  data  is  given  which  is 
harmonic  in  space  (i.e.,  of  the  form  e     )  and  that  the 
resulting  solutions  are  of  the  form  a(t)e     .   The  vector 
a(t)  will  be  shown  to  be  unique.   The  method  of  proof  will 
be  seen  to  be  applicable  to  solutions  whose  Fourier  trans- 
forms exist  for  some  time  T  >  0.   In  this  case,  the  time 
evolution  of  the  Fourier  transform  would  be  shown  to  be 
unique . 

We  write 


u(t  ,xj  =  u(t )e 


B(t,x)  =   B(t)e^^''' 

f  (t,x;  [x,v)    =  f  (t;  ii,v)e^^'^ 


and  substitute  in  eq.  (ll)-(l3)  to  obtain 


+ 

+  r  +  +  f 


(16)       0  =  (^  +  lk^v)f~  +  lk^B^U[ri]a"  -i^)^'^^^-   v  ^ 


+ 


f         t   t      2\ 
+  I(k2r2  +  k  r  )  g°  +  Ik^a"  f"   <v  > 

o 


SB. 
(17)      0  =  ^  +  l(k2r'2  +  k^r^) 


°  -^^2,3  "  ^V2,3 


b2 


(18)      0  =  ^  r^^^  -  lk,v%^^  +  Ik^^^  c\   +  ^  ik2^3<^i> 

where  we  have  used  k  =  (k  ^.kp  ,k^),  and  the  quantities  <(  ix)i 
and  \v  y  appearing  here  are  the  same  as  those  of  equations 
(15)   modulo  the  factor   e   "  '.     If  we  assume  two  deriva- 
tives for  the  dependent  variables,  then  we  may  eliminate  r^ 
and  r  to  obtain,  for  B-,  ,  the  equation 

2  2    /,  2    ,  2x  2K„        o  /,  2    ,  2 


(19)       (^  +  ^k^v^  +  {14   +    k^)c^  )B.  .  --°  (k^  +  <)<^) 
St    I  ^  ^o 

Equation  (19)  can  also  be  obtained  by  Integrating  equations 
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(18)  with  respect  to  time  and  using  equations  (l6)  and  (17) 
to  prove  the  existence  of  the  second  derivative.   Solving 
equation  (19)  for  B-^  in  terms  of  (m-),  we  have 

i(ot)      -±{\t) 
(20)      B^  =  b^e    ^  +  b^e 


dT  sin[fi^(t-T)](^i(T)) 


^o      k 

o 

where  we  have  written 


ll/2 
(21)      0^  =  ^v^k^  +  (k^  +  k|)c2 


and  b-,  and  bp  are  homogeneous  of  degree  one  in  B-j^(O)  and 
^B-j/St  ,_^,    the  latter  being  given  by  the  initial  data  of 

rp(0)  and  r^(0).   In  a  later  section,  it  will  be  shown  that 

2 
if  V  <  0,  the  initial  value  problem  is  not  well  set  for 

2 
the  linear  system.   Hence,  we  consider  only  the  case  v  >  0, 

_\ 
which  guarantees  that  Q,  is  real  for  real  k.   It  should  be 

noted  that  if  Q,     is  not  real,  the  proof  of  uniqueness  is 

still  valid  with  minor  modifications  in  the  estimates  which 

will  be  made.   However,  the  estimates  of  grov;th  in  time  are 

very  different.   The  reality  of  Q,    will  be  seen  to  yield 

linear  exponential  bounds  on  the  functions,  whereas  the  well- 

2 
posedness  theorem  for  v  <  0  vjill  be  shown  to  imply  otherwise. 

That  is,  in  this  case,  the  problem  cannot  be  solved  by 

Laplace  transforms  in  the  space  of  solutions  having  Fourier 
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transforms  in  space. 

Returning  to  the  proof,  we  replace  the  expression 
l(kpr  +  k,r  )  in  eq.  (l6)  by  -BB-|/St  from  eq.  (17),  thereby 


I' 


+ 


obtaining  equations  for  f~  in  terms  of  B-,  .   But  B-,  is  known 
as  a  functional  of  f~,  and,  therefore,  we  obtain  from  eq.  (l6) 
expressions  of  the  form 


(21)       (^+  ik^v)f-=G-(f^(T),f-(T)) 

where  G(f  ,f~)  is  a  linear  functional  of  the  vector  (f"*',f~) 
and  is  composed  of  integrals  thereof.  We  write  eq.  (21)  in 
the  equivalent  form 


(22) 


^  4.       -ik-jVt 

f-(t;kL,v)  -  f-(0;ki,v)e 


.t 


+ 


•ik  V(t-T)  + 

■        1         G-(f+(T),f-(T))dT 


and  have,  upon  expanding  G  under  the  integral  sign,  the 
lengthy  expression 

t 


+ 


+    -ik-,vt 


(25)      f-(t)  =f^(0)e 


1 


+ 


+  +       ,  p   .  -ik  v(t-T) 


o 


+ 


dT  e 


ik-,  V  (t-x) 


^1 


+ 
f  ~      ifi,  T      -ifi  T 

[ifi,  ^  (b,e   ^   -  b.e   ^  )] 


k  B   ^  1 
o 


'2 


+ 


r  -ik_v(t-T) 


+ 


dr  e 


+ 


J 
o 


[-ik  ',  ([nJa--  li)f„  ^^  -V  w^  (b.e 
1^  o,v     B^J   1 


in,  T 

k 


-ifi,  T 
+  b.e   ^  )] 
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(23)  ^ 

f^        -Ik  v(t-T)     f     +     +      f:  /B^  (k^  +  k!) 


B   p     0., 
o  ^o     k 


X 


d^  sin[a  (t-^)]<h(^)>] 


+ 


-ik-,v(t-T)      f^        o     o    r 

dT  e   1      [_^B(k2+k?)    d^  cos[fi,.(T-^)]<n(^)>] 


'k 


V   S    // 


We  now  reproduce  the  quantities  \  M.(t  )>  and  <"  v  (t)\  from  eq. 

(23)  by  first  multiplying  by  |i.  the  equations  for  f  and  f 

2     +        2     — 
and  adding,  and  then  multiplying  by  v  (l/p  )  and  v  (l/p  ) 

respectively  and  subtracting.   Integrations  with  respect  to 

M-  and  V  then  yield  expressions  of  the  type 


(21+) 


\M-(t)>l    <^  P^^^t)   +  a^    I    dC  |<v2(T)>|dT 

o  o 


l^ 


,t       X 


+   b. 


V- 


dC 

o  o 

t         r  T 


<^^(T)>|d' 


r 


+     C 


IJ- 


dC  dT      I       I  ([1(1))  I  d^ 


O  O 


J 


(25)  Kv2(t)>|  <p^^^)(t;  + 


dC 


V      J  J 

O  O 


<v2(T))|dT 


+  b 


v2. 


-t         A 


dT 


h>(^)>|de 


'o  o 

t  C  T 


+   c 


v2  J 


dC        dT    j     |/n(|))|d^ 
000 
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In  obtaining  (2^)  and  (25)  the  exponentials  In  (23)  were 

majorized  In  absolute  value  by  unity.   The  quantities  a  , 

b  ,  c  ,  a  p   b  p,  and  c  p  are  maximal  estimates  for  the 

^        ^        V  '   V   _(.     v"^  _^ 
Integrals  of    f  and   f~   with  respect  to  the  moments 

which  occur  together  with  the  coefficients  of  equation  (23) 

which  multiply  them.   The  expressions  P^   (t)  and  PA,"^   (t) 

are  second -degree  polynomials  which  are  positive  and  are 

homogeneous  of  degree  one  In  f~(0;M-,v),  b-,  and  bp . 

We  define  the  quantity  m  =  max  |  a  ,b  ,c  ,a  p,b  p,c  p} 

^     ^      ^      ^      V         Y         V    ' 
and  the  functions  X(t)  and  Y(t)  as 


X(t) 


f 


t   X        A 


dT 


l<^^(^)>|dl 


(26) 


■J  J  ^ 

000 


Y(t)  = 


dC 

dC   dT  |/v^(t)\ 

J    J 

o    o 


together  with  the  functions  R(t)  and  S(t)  as 


R(t) 


^\ 


(27) 


at- 


I  -  m  Y(t)  +  (^  +  l)X(t) 


2 

S(t)  =  ^  -  m[Y(t)  +  (^+  l)X(t) 
St      I 


Then,  in  view  of  the  definitions  of  X(t),  Y(t),  and  m, 
reference  to  equations  (24)  and  (25)  yields  that 


(28) 


R(t)|  <  P^^^t) 
S(t)|  <  P^^^)(t)  . 
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We  now  solve  equations  (27)  and  use  (28)  to  obtain  estimates 

for  the  functions  X(t),  Y(t)  and  their  derivatives.   These, 

o 
In  turn,  will  yield  bounds  on  KkL(t))>|  and  |<^v  (t)^|  via 

Inequalities  (24)  and  (25).   We  reduce  (27)  to  a  first-order 

system  by  defining  the  vector 


(29) 


p(t)  =  1|,  M,  X,  ai,  y' 

\St   St     ht        J 


and  obtain 


(30) 


^ 


P(t) 


r 


m 


0 

1 

1 

0 

1 

r       \ 

R(t) 

1 

m 

0 

0 

0 

0 

0 

0 

1 

m 

0 

0 

0 

p(t)  = 

0 

0 

1 

1 

0 

1 

s(t) 

0 

0 

0 

1 

m 

0 

.°  . 

together  with  the  Initial  data 


(31) 


P(0)  =  0 


From  eq.  (30)  and  (31)  we  easily  obtain  that 


(32) 


^      5 


Pj_(t)  <  (3m+l) 


r 


1=1 


P^(T)dT  + 


R(t)  +  S(t)  dT 


^o 


1  =  (1,2,3,4,5) 
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or 


(33) 


r     3 


P, (t)  <  51  F. (t)  <  5(3m+l)   y  F  (T)dT 
1  1  i=l   ^ 


+  5 


(R(t)  +  S(T))dT, 


o 


Defining 


(3^)      H(t)  = 


51  Pi('^)dT 


1=1 


we  have 


(35) 


§1  -  5(3m+l)H(t)  <  5 


H(0)  =  0 


(T(t)  +  S(T))dT 


o 


which  Implies 


(36)      H(t)  <  5 


^5(3m+l)(t-T) 


J 
o 


R(^)  +  S(e)   dl 


o 


Substituting  In  equation  (32)  and  using  the  bounds  from 
(28)  we  finally  obtain 


(37)      Pi(t)  <  5(3m+l) 


^5(3m+l)(t-T) 


P^^)(e) 


+  P^""  ^O  d^ 


since  P^^^(^)  and  P^'^  ''(4)  depend  linearly  and  homogeneously 
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on  the  Initial  data,  this  yields  the  uniqueness  of  the  Pj_(t), 

o 

and,  consequently,  of  |(l^(t)>!  and  |<v  (t)>|  via  eq.  (24) 
and  (25).   The  polynomial  bound  on  P^^''  and  P^^  ^    in  turn 
implies  a  linear  exponential  bound  on  the  F^  and,  as  a 
result,  on  KM-(t)^|  and  |(v  (t)>|.   Using  this  bound  on 
|<;^|i.(t))>|  ,  equation  (20)  yields  a  similar  bound  on  B-j^(t) 
together  with  uniqueness.   If  we  integrate  equations  (17) 
and  (18)  and  eliminate  r^    -2^  we  obtain  the  necessary  esti- 
mates for  B^  ^.   Eliminating  B^  ^,  we  obtain  the  estimates 

for  r^  v- 

^'^  +      ± 

To  recapitulate,  we  have  shown  that  if  f  and  f 

decay  sufficiently  fast  so  that  the  moments  which  appear  are 

absolutely  integrable,  then  solutions  which  are  harmonic  in 

space,  or  which  have  Fourier  transforms  for  some  time  T  >  0, 

are  unique  within  the  class  of  such  functions.   Moreover,  if 

-^  {BJii     +  P°  -  P°)  >  0,  then  the  time  rate  of  growth  of 
Pq      o 

these  solutions,  on  the  transforms,  is  bounded  by  a  linear 
exponential  function  of  time  which  is  a  function  of  the  wave 
vector  or  of  the  transform  vector. 
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TV.   Stability 

A.   The  Transformed  Equations  and  a  Well-posedness  Theorem. 

Before  proceeding  to  an  analysis  of  wave  propagation 
and  stability,  we  will  prove  a  well-posedness  theorem  for 
the  Initial  value  problem  associated  with  the  system  of 
equations  (11)-(13). 

We  say  that  the  problem  is  well  posed  if  and  only  if  a 
unique  solution  exists  for  the  initial  data  and,  moreover, 
the  solution  depends  continuously  on  the  data.   That  is, 
if  s(t;l(x))  is  the  solution  corresponding  to  the  initial 
data  I(x),  then  we  require  that  for  every  e  >  0  there  should 
exist  a  t^  >  0  and  a  5  >  0,  such  that  if  |f(x)|  <  6  and 
t  <  t^,  then  I s(t;l(x)) I  <  e. 

The  theorem  we  prove  is  the  following:   if  initial 
data  is  given  which  is  harmonic  in  space  with  wave  vector 
k,  and  if,  corresponding  to  k,  there  are  solutions  of  the 
form  e"^^'^^^)^-^^*^^)  with  Im(a3(k))  >  0,  then  the  initial 
value  problem  for  our  system  is  not  well  posed. 

Suppose  that  such  co(k)  existed,  then  substitution  of 
the  solution  into  the  system  yields  an  algebraic  system 
which  is  homogeneous  of  degree  one  in  the  vector  k=  (co(k),k). 
Thus,  if  we  replace  k  by  A  k  and  Go(k)  by  A^a)(k),  then  the 
equations  are  still  satisfied.   Now  suppose  that  an  eigen- 
vector of  the  system  associated  with  a3(k)  is  the  Initial 
data  amplitude  vector  v^,  some  of  whose  components  of  course 
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contain  the  parameters  \x   and  v.   Then  If  we  associate  with 

A  k  the  Initial  data  vector  v  /A  and  the  frequency  A  cu(k) 

2 
and  consider  the  limit  as  A  — >  +00,  we  obtain  Initial  data 

which  Is  arbitrarily  small,  but  solutions  which  are  arbi- 
trarily large  for  any  finite  time.   This  completes  the  proof. 

In  particular.  If  we  take  the  plane  of  the  vector  k 
to  be  the  x-y  plane,  and  the  Initial  data  vector  such  that 
all  dependent  variables  are  zero  except  r^  and  B^,  then 
the  system  becomes 


(38) 


0  =  -vT^  -  Ik-,  r^  =^  0  =  cjoB^  -  k^  r., 
ot      1  p  ^    11 


0  =  ^^  -  Ik,  v^B^  =^  0  =  a3r^  -  k-,  v^B^  , 
ot      1        p      ^  513 


so  that  If  k-,  /=   0,  co(k)  =  t  k  v  sufflcles  for  the  existence 

of  non-trlvlal  solutions.   Therefore,  the  posltlvlty  of 

2 
V   IS  necessary  for  well  posedness  and  exponential  stability. 

2 
In  passing,  we  note  that  this  Inequality  (v   >  O)  was 

obtained  by  Liist   as  a  condition  for  stability  from  a 

macroscopic  plasma  model. 

As  a  result  of  the  above  conclusions,  we  will  assume 

2 
in  the  following  that  v  >  0,  and  use  the  results  of  the 

uniqueness  theorem  to  solve  the  system  of  equations  (I6)- 

(18)  by  Laplace  transform  analysis.   If  p  Is  the  transform 

variable,  we  obtain 
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(38) 


(^+  v)?-  +  S^[([^]a- 


P 


+ 


^i)f    -  V  „ 

o,v     B 


!i) 


o 


k. 


k- 


^P  ^    '^^  '^    o    "'"  - 

+  (k7^2  +  kf^j)  B-+«'^o,v-/  -   ik. 
1       1       o  i 


tf^  <-2\  ^  fiLo) 


(39) 


P 


B. 


Ik^  ^2,3 


2,3 


Ik. 


(40) 


P 


ik  ^2,3 


2^ 
V  B2^3  + 


k 


2^2^ 
r^  c'^B.  + 

k-,      1 


Bo  '^2,3 


o 


k- 


<^) 


''2.3'°: 

Ik, 


where  we  denote 


(41) 


Q(p)  = 


-pt 


,(t)dt 


J 
o 


Ik  •  X 


if  we  are  referring  to  solutions  of  the  form  Q(t)e     ,  or 


(42) 


;(p,k)  =   e 

4) 


-pt  ~ 


(t,k)dt 


if  Q(t,k)  is  the  Fourier  transform  of  Q(t;x). 

In  our  discussion  of  stability,  it  will  suffice  to 
keep  the  vector  k  fixed.   If  it  is  a  Fourier  transform 
variable,  then  all  remarks  pertaining  to  the  "solution"  of 
our  discussion  will  be  applicable  to  the  "Fourier  transforms" 

Thus,  we  now  take  the  plane  of  k  to  be  the  x-y  plane. 
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The  wave  vector  k  is  assumed  to  make  an  angle  0,    not 
equal  to  t  7r/2  with  the  vector  B ^ .   (The  special  case 
6   =   7r/2  will  be  considered  separately  and  will  be  solved 
exactly.)   After  defining  the  variables  z  =  -p/lk-^   and 
■X  =  tan  0,   we  have 

(^3)       (v-z)f-  +  B^[([n]a--lx)f^^J(v-z)B^ 

--fo,v<->=il^(^'(°^^'-)  -^'o) 

B  (0) 
-zB-,  +  Ar. 


(44) 


a     2    ik^ 

B    (0) 


'2,^  2,3     Ik^ 


b!  ...   r^(o) 


-zr 


v^Bo  +  ^c  B,  +  —  A<n>  =  -^ 


(^5) 


2    '  ^2      ^1    p    '"^'     ik- 

^o 


2-     ^%(0) 
-zr-^  -  V  B-  - 


3      3    ik- 


Dividlng  eq.  (43)  by  (v-z)  and  multiplying  by  the  appro- 

2  ± 
priate  kernels  (i.e.,  v  /p   or  |J-)  and  combining  to 

reproduce  the  quantities  <^l^)>  and  <^v  ^,  we  obtain  the 

following  system  for  the  vector  s(z)  =  (vM-)>,s,v  /•,B-,,Bp,rp 


(46)       A(z;A)a(z)  =  ^t{0)    , 
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where  j(0)  Is  a  vector  given  by 


[J- 


g^(Q)+s"(Q)  d^dv 


^ 


v-z 


V 


S^(0)/p':  -s'(0)/p, 


V-Z 


-dio-dv 


(^7)    j(0)  = 


6^(0) 
62(0) 
r.(0) 


6^(0) 


and  A(z;A)  Is  a  7  x  7  matrix: 

(^8)    1   -Vl-|fo,v)A-((^^V(^o,v)x 


A  = 


1 


0 


T 


[[.]T-z%(^X^l 


0 

0 

0 

0 

< 

0 

-z 


0 


0 


0 


Ac' 

0 
0 


;  g-(o)  =  f-(o) 

B  (0)  + 
— f  " 

B.     o 


"^ 


0    0    0    0 


0    0    0    0 


0 

A 

0 

0 

z 

-1 

0 

0 

2 

V 

-z 

0 

0 

0    0   -z   -1 


0   0   -V   -z 


J 
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in  which  the  following  notation  has  been  employed 


I   ^  (h) 


00      -\-co 

r       ^ 

dM- 

J    J 

O      -00 


k 


^'h(^i,v)  ^^ 


v-z 


(48-)      q^  =  q+  -  q- 


q^  =  q"*"  +  q" 


T  =  (l+|al^)  -  z   I     (Ia|  f^yp^)^ 


A  calculation  yields  that 


(49)       det  A(z;A) 


-  z(z^-v^)D(z;a2) 


with  D(z;A  )  given  by 

pO        g2 

(50)       D(z;a2)  ^  t(z2-[v2+a2(c2+  ^)  ]  -  A^  ^  I  ^(fQ  ^)v) 

^O  ^O    l-L      ' 


-zV^(lalp)-   l^(f^yp^)^' 


We  note  that   A  may  be  v;ritten  as 


/" 


(i^8')     A  = 


M 


0 


^ 

0 

-z 

-1 

2 

-V 

-z 

with  M  a  5  X  5  matrix  with  determinant 
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(51)      det  M  =  -  zD(z;?s^) 


If  we  assume  the  Hblder  continuity  and  absolute  Integra- 

°°     If  "'" 
blllty  of  /  M-  f"   di-L  (k  =  0,1,2),  then  all  the  elements 
o     o,  V 

of  A(z;A)  are  analytic  In  z  in  a  half  plane  and  continuous 

15  l4 
up  to  the  real  line.  -^ '  We  note  that,  in  general,  the 

elements  cannot  be  continued  analytically  across  this  line, 
and,  therefore,  if  k  is  a  transform  variable  it  must  be 
restricted  to  a  half  space,  k^  <  0  or  k.^^  >  0.   The  qualita- 
tive results  for  stability  are  the  same  in  both  cases, 
since  if  p  yields  an  instability  for  k-j^  >  0,  -p  yields  the 
same  result  for  -k-,  ,  as  an  inspection  of  D(z;A  )  readily 

reveals.   Thus,  we  take,  without  loss  of  generality,  k-j^  >  0 

2 
and  our  criterion  for  stability  becomes,  assuming  v  >  0, 

2 
that  D(z;A  )  shall  have  no  zeros  in  the  upper  half  of  the 

z-plane,  which  is  its  domain  of  analyticity  (assuming  the 
absolute  integrability  of  the  quantities  mentioned  above). 
It  actually  suffices  that  they  belong  to  L  ,  p  >  1. 

Before  proceeding  to  the  investigation  of  stability 
we  note  that   had  we  sought  plane  wave  solutions  for  a 
fixed  wave  vector,  the  complex  frequencies,  oo,  of  propaga- 
tion would  have  appeared  as  roots  of  the  determinant  of  the 
matrix  A,  with  z  =  t  co/k.,  .   However,  it  would  have  been 
Impossible  to  simply  discover  real  roots  of  D(z;A  )  since 
these  would  give  rise  to  improper  integrals  (i.e.,  v-z 
would  vanish  for  some  v  in  the  range  of  integration)  and  a 
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special  Interpretation  of  such  Integrals  would  have  been 
necessary."''^  An  evaluation  of  the  Inverse  Laplace  transform 
by  Integration  along  a  path  in  the  upper  half  plane  avoids 
this  problem,  and  in  yielding  the  complete  solution  to  the 
problem  uncovers  all  such  modes  of  propagation. 


B.   Stability  Criteria 

In  the  following  derivation  of  stability  criteria  for 

9 
the  steady  state  configuration,  we  follow  0.  Penrose   and 

exploit  the  principle  of  the  argument  for  analytic  func- 
tions.    The  statement  of  the  theorem  we  employ  is  the 
following:   If  f(z)  is  analytic  in  a  simply  connected  domain 
^  and  continuous  up  to  its  boundary  S^,  then  the  number  of 
zeros  of  f(z)  in  ^  is  equal  to  the  index  of  the  image  of  ^3' 
with  respect  to  zero  when  hp   is  taken  v;ith  u    to  its  left  in 
the  z -plane. 

With  a  view  toward  the  application  of  the  theorem,  we 
consider  functions  of  the  form 


(52)      f(z)  =  z^a(z)  +  b(z) 


with  a(z)  and  b(z)  analytic  and  bounded  in  the  upper  half 
plane  and  continuous  up  to  the  real  like.   We  further  require 
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that  for  z  belonging  to  this  half  plane 


11m   a(z)  =  a  =  real  positive  number 
I  z  I  — >  «> 
(53) 

11m   t)(z)  =  b  =  real  number. 

(z|->  00 


Also,  for  z  real 


a(z)  =  a(-z) 

(5^)  

b(z)  ^   b(-z) 

where  the  i   denotes  the  l-conjugate. 

We  now  consider  the  mapping  of  the  domain  bounded  by 

1* 
the  curves  C-,:z  =  Re   ;  0<*<7r  and  02  =  2  =  x  +  lO; 

-R  <  X  <  +R  with  R  »  1  (see  Fig.  l).   As  a  result  of  eq. 

(53),  C-,  Is  mapped  Into  a  large  arc,  C-,  ,  which  approximates 

a  circle  crossing  the  real  axis  In  the  Image  plane  near  -co. 

As  <i>   moves  through  Increasing  values,  the  direction  taken 

by  C-|  Is  counter-clockwise  (see  Figs.  2,3,4-).   As  a 

result  of  the  conditions  Imposed  on  a(z)  and  b(z)  In  eq. 

(5^)  we  are  assured  that  the  Image,  Cp ,  of  C^  crosses  the 

real  axis  when  z  =  0,  and  the  Image  for  z  negative  Is  the 

conjugate  of  that  for  z  positive  (i.e.,  they  are  reflections 

of  each  other  In  the  real  axis  of  the  Image  plane.   However, 

the  direction  vectors  of  the  tangents  are  Images  of  each 

other  In  the  Imaginary  axis  of  the  plane.) 
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Pigs.  1-4  are  chosen  to  Illustrate  the  nature  of  the 
possibilities  that  may  arise  and  to  motivate  the  derivation 
of  the  stability  criteria. 

Fig.  1  simply  Illustrates  the  domain  O"  (striped  area) 
and  locates  some  points  which  are  symmetrically  placed  on 
the  real  axis,  so  that  the  direction  taken  by  the  mapping 
may  be  followed  easily. 

In  Pigs.  2-4,  the  vertical  axes  are  not  drawn  so  that 
the  nature  of  the  mappings  may  be  distinguished  when  f(0)  >  0 
and  f(0)  <  0. 

In  Fig.  2  we  have  the  simplest  mapping.   It  is  one-to- 
one  with  the  image  of  the  real  axis  crossing  the  real  axis 
in  the  image  domain  only  at  z  =  0.   We  note  that  if  f(0)  <  0 
the  index  of  the  image  of  SD  with  respect  to  zero  is  zero, 
and  if  f(0)  >  0  the  index  is  one.   Hence,  in  this  case,  if 
f(0)  <  0  then  f(z)  does  not  take  on  the  value  zero  in  the 
upper  half  plane,  but  if  f(0)  >  0  it  assumes  this  value  once. 

Fig.  5  illustrates  a  slightly  more  complicated  mapping. 
Here  the  image  of  Cp  crosses  the  real  axis  at  z  =  "t  x.,  in 
addition  to  z  =  0  to  yield  a  mapping  which  is  not  one-to-one. 
We  note  that  if  f (O)  >  0  then  the  index  of  the  boundary  image 
with  respect  to  the  origin  is  one,  and  f(z)  assumes  the  value 
zero  just  once  in  the  half-plane  Im  x  >  0.   If  f (O)  <  0  and 


f (-  X. )  >  0  then  the  index  with  respect  to  the  origin  is 
and  f(z)  assumes  the  value  zero  twice  in  this  half  plane 


two 


When  f (O)  <  0  and  f (t  x, )  <  0  we  have  an  index  of  zero  with 
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respect  to  the  origin,  so  that  f(z)  Is  never  equal  to  zero 
in  the  upper  half  plane. 

The  Inference  we  are  led  to  by  the  above  considerations 
Is  that  when  f(0)  <  0,  f(z)  will  take  the  value  zero  2n 
times  (n  =  0,1,2,...),  but  when  f(o)  >  0  it  will  assume  this 
value  2n+l  times  (for  the  same  range  of  n).   We  now  prove 
this  theorem.   Since  f(t  oo)  are  real  and  positive  and  f(0) 
is  real,  the  Images  of  the  half  lines  x + 10  >  0  and 
X  +  10  <  0  induce  an  index  with  respect  to  the  origin  of  -^n 
(n=-l,0,l,2, . . . ) .  If  f(0)  <  0  then  n  is  odd.   If  f(0)  >  0 
then  n  is  even.   As  a  result  of  the  symmetry  of  the  images 
of  these  two  half  lines,  these  indices  are  the  same  and 
supplementary.   That  is,  if  f (O)  <  0  then  the  index  due  to 
the  real  line  is  2  x  -^n  with  n  odd,  while  if  f(0)  >  0  then 
the  index  is  2  x  ^n  with  n  even.   Thus,  adding  the  index  one 
due  to  the  image  of  C-,  ,  we  have  proven  the  theorem. 

From  the  arguments  preceding  the  theorem,  one  might  be 
tempted  to  conclude  that  any  additional  crossings  of  the 
real  axis  in  the  image  plane  must  also  occur  at  negative 
values.   However,  as  is  shown  in  Fig.  ^,    this  is  not  the 
case.   Here  we  see  that  if  f(0)  <  0  together  with  both 
f(t  X-,)  >  0  and  f(t  Xp )  >  0,  then  the  index  of  the  image  of 
SD  with  respect  to  the  origin  is  still  zero,  so  that  f(z) 
has  no  zeros  in  D.   It  is  obvious  that  such  a  condition 
together  with  f(0)  <  0  would,  however,  be  sufficient  for  f (z) 
not  to  have  zeros  in  the  half  plane. 
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2 
We  now  return  to  the  function  D(z;A  )  and  Impose  the 

+ 
conditions  that  f  {\i,v)    should  be  even  in  v,  and,  in  partlcu- 

2 
lar,  that  they  should  be  functions  of  v  .   Then  we  have, 

together  with  the  conditions  imposed  on  a(z)  and  b(z)  above, 

that  the  remarks  made  pertaining  to  zeros  of  f(z)  above  are 

applicable  to  D(z;A  )  (see  Appendix).   However,  exponential 

instability  is  equivalent  to  the  existence  of  zeros  of 

2  2 

D(z;A  )  in  the  upper  half  plane  of  z,  assuming  v  >  0,  and, 

hence,  for  exponential  stability  of  the  plasma,  we  must  have 


(55) 


o    2 

0  >  D(0;A-)  =  -(1+  \a\^)\v^+-h^ic^+~   +  ^  P-V. 

"I  Po   Po 


X 


\^  ^!^a^a.) 


where  P.V.  denotes  principle  value.   With  the  hypothesis  f- 

2 

are  functions  of  v   rather  than  v;  the  principle  value  is  not 

necessary  since 


(56) 


P.V, 


M- 


.-  (^o,v)v  dixdv  =  2 


H^(f 


)d|i.dv 


o,v 


Since    v     >   0  and  we   want  D(0;A^)   <   0  for  all   finite   A^ ,    eq. 
(55)   becomes 


(57) 


c      +  —  +   -^  2 
Pq        Po 


r 

f 

dn 

« 

^ 

.+  OJ 


dv   ^i^(f        2)      >   0 
o,v     2 


1 


B 


o 


o. 


p       (—  +Px) 
^o        o 
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Eq-  (57)  presents  us  with  a  relatively  simple  first  test 
for  the  stability  of  the  plasma.   In  particular,  it  is 
simple  to  show  that  all  spherically  symmetric  functions 
(i.e.,  functions  of  the  form  f~  -  fo^l-^^o  +  ^v  )  )  satisfy  the 
inequality.   For  such  functions  we  have 


OjV        o     '^ 


Substituting  eq.(58)  into  (57)  and  integrating  by  parts 
with  respect  to  1-L  yields 


(59)       ^'-r-f^o, 


O  ^00       +00  o        ^2 

^d^x    (f^)^dv  =  c2-^  =  ^^>  0, 

J  "^o  "^  o   o 

O       -00 


We  note  that  for  such  functions  p°~  =  p°  ,  so  that 
V  =  -—  B  /m-  and  the  two  necessary  conditions  for 
stability  degenerate  to  a  single  one. 

While  eq.  (57)  presents  us  with  a  simply  verifiable 
necessary  condition  for  stability,  something  more  than 
this  may  be  obtained  from  it  by  introducing  dimensionless 
variables.   Thus,  we  write 


(60)      fo(kt.v2)  =  a+  g^(b^n,  c^v^) 

and  evaluate  a^,  b+ ,  and  c_^  in  terms  of  p~,  p°~,  and  p°~. 
Doing  so  we  obtain 


See  appendix. 
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P„    //  Soly^x  jdydx 


(61) 


t  +2 

P^        //yg^(y.x    )dydx 


x  o      o± 


2 


Pi"      //   &Jy'^    )dydx 


a+   - 


Bjp^)5/^      (//x\^(y,x2)dydx)2//yg^(y,x2)dydx 


Ot  /       0+  N  2 

Px     (P„     ) 


(//   g3(y,x2)dydx)5/2 


2    2  2 

Substituting  y  =  b^ii.;    c^v     =  x      In  eq.    (57)    we   obtain 


(62) 


2    Sf 


1)^0.^0 


0<:^^^^+2p,  +    2(     ^^ 


,    o^2      /x    g  dydx   /   y     ^  dydx' 


(/y  gdydx)' 


2 


If  we  now  assume  that  Pj_   =  2  p_i_  and  p  "  =  i^  p„  together  with 
~~(y^x  )  =  h(y,x  ),  we  obtain,  using  the  Identity  2x  — p  =  -^ 


by. 


and  integrating  by  parts 


B 


o2 
1    Px 


(63)  0  <    (-^+    2p^)    -  ^-^Ti 


'\^ 


o 


3      o 

P 
II 


with   T)   given   by 


(6^) 


n  = 


(/   x"^  ^dxdy)(/   y2   ^dxdy) 

ox ox 

f  r      2      ah      ^    ^    n2 
(J    X   y  — 2   dxdy) 

Bx 


,10 


Liist,      using  a  macroscopic   model,    has  derived   eq.     (63)   v;lth 
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2  2 

Tj   =   1 .      If  h(y,x    )    is  monotone   down   in  x     for  every  value 

2  I  2  1  ~  2  \  ^ 

of  y,    then,    taking  h-,    =  x    |^h/Sx    |^   and   h^   =   y|Sh/Sx    |^   and 

using  Schwartz's  inequality,  we  have  r;  >  l  so  that  this 

class  of  functions,  at  any  rate,  for  stability  the  guiding 

center  theory  imposes  greater  restrictions  on  the  pressures 

than  the  fluid  theory  does.   Agreement  is  obtained  in 

limiting  cases  of  functions  belonging  to  this  class.   In 

particular,  if  we  take  the  sequence  of  functions 

h.  =  e-J^(y-yo)' — L^       y  >o 

J  1+  (x^)J       ° 

and  define  t]  .  as  the  value  of  r\   when  h.  replaces  h  in  (64), 

J  J 

then  we  obtain 


lim  T)  .  =  1 
j  — >°°   "^ 

If  we  are  not  restricted  to  monotonically  decreasing 

2 
functions  in  x  ,  the  quantity  rj  may  be  less  than  one  or 

even  negative  (in  which  case  eq.  (65)  is  always  satisfied), 

as  the  following  example  shows.   We  choose  h(y,x  )  such  that 


Then 


2  x^ 

h(y,x^)  =  H(y)  ^ 

(1+x  )"^ 


ax^        i  (l+x^)"'+^ 


-  31  - 


Referring  to  eq.  (62),  since  H(y)  >  0  [h{y,x    )  is  an  expecta- 
tion function],  it  suffices  to  show  that 


(65) 


(l+x2)-+l   ^"> 


o 


Letting  x  =   tan   0   in  eq.    (65),   we    obtain 
7r/2 


m  cos^^'e    -    (m-l)cos^^"^"^^0      dS 


TT     1-3-  •  •(2m-3)    ( 


IT  1  2.4.  .  .  (2m-2 


>    0 


In  Pigs.  5  and  6  we  indicate  the  ranges  of  p   and  p_i_  which 
are  required  for  stability  by  eq.  (63)  and  the  condition 
V   >  0  when  t]  >  0 .   We  note  that  if  rj  >  6,  then  for  a 
fixed  magnetic  field  strength,  B  ,  there  exist  maximal 
pressures  Pj_  and  p   such  that  if  Pj_  >  p_i_  or  p^^  >  p^^  , 
stability  is  impossible.   On  the  basis  of  such  a  predic- 
tion, it  might  be  anticipated  that  for  all  the  functions 
which  we  could  admit  to  consideration  (i.e.,  positive  and 
vanishing  at  infinity  in  the  appropriate  |x-v  half -plane) 
r\   must  be  less  than  six.   The  following  counter-example 
shows  otherwise: 


h(y,x2)  =   1 


-X 


(1+H) 


621z2  >  6 
m-5 


m 


if 


m  >  4 


-32- 


Figures  7  and  8  are  plots  of  the  temperature  ratios  against 

P  =Px/(B>o)- 

Before  taking  leave  of  this  criterion  we  note  that 

all  spherically  symmetric,  and  consequently,  all  two- 
temperature  Maxwelllan  distribution  functions,  have  the 
value  six  for  tj  . 

As  remarked  above  during  the  discussion  of  zeros  of 
f(z).  It  Is  not  sufficient  to  Insure  exponential  stability 
that  the  two  conditions  derived  thus  far  should  be  satisfied 

That  Is,  D(0;A  )  may  be  negative,  and  there  may  still  exist 

2 
zeros  of  D(z;A  )  In  the  upper  half  plane.   Therefore,  we 

now  turn  our  attention  to  the  derivation  of  a  condition 

which,  together  with  the  two  necessary  conditions  for 

stability,  will  guarantee  the  stability  of  the  equilibrium 

configuration.   This  condition,  not  surprisingly,  will 

depend  on  microscopic  properties  of  the  plasma. 

2 
The  sufficient  condition  we  impose  is  that  D(z;A  ) 

should  never  be  real  and  non-negative  for  z  real.   As  noted 

2 
previously,  together  with  D(0;A  )  >  0,  this  implies  that 

2 
D(z;A  )  is  never  zero  in  the  half  plane. 

2 
Accordingly,  suppose  that  D(z;A  )  is  real  for  some  z 

which  is  real  and  different  than  zero.   Then  for  such  a 

2 
value  of  z,  assuming  A  /^  0  and 


/  0  , 


^l(—  f    ).d\± 


v=z 
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we  have,    from  Im  D(z;A    )    -   0,    that 


(66)  i„(7^f.  J      -I 
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"^^  Po      °'^    A 


A     — ^      (ap    )  ^  P^      o,v^A 

^  o 


X    \{Re   T(z))(-A 


•r2 

2    ^o 


o 


^  (f.  .)  dn) 


o,v'^ 


o 


-   z 


a 


o  2 

2      2      2  /    2      -P-"- 


p         o,vj  ^  ^  p  p  2^o,vV 

^0/£  "^O^OiJ,  '       ^ 


where  x  ,  =  P.V.  I  ,  for  z  real.   Substituting  this  into  the 


^ 


M- 


expresslon  for  Re  D(z;A  )  when  z  is  real,  we  obtain 


o        ^2 
(67)  Re   D(z;a2)    =    (Re   T(z)    z2-v2_a2(c2+^+  ^i   ^(f         )     . 


B 


2 


/  00 


^2^2    ^o      2, 
•^o 


a 


I  f .   ,\ 


o,v 


I   o 


V 


aix 


+ 


IJ-    (f        )   di-L-  ap 
^    o,v^  '    ^o 


i-ri'^ 


-    (Re    T)A 


2       o 


2, 


o        ^2 
2/2       2^2/2       P-L 


o 


(^i   f         )    dH  +    z^    z^-v^-A^    c^+  —  +   -^  f  ^(f^      ,, 
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a 


U 


(^    f  )     d^LJ 

Po       ^'^   Z         ' 


o 


^^'^'S"-o)^-l      [[{-ffl 


2 


dp. 


-   3^    - 


If  the  first  product  appearing  In  eq.  (67)  Is  negative, 
then  It  suffices  to  have  at  v  =  z 


(68) 


I  -  ^ 


2 


00 


o  /A 


'V^Vl- 


o 


^  (f^   )   d^i 


We  now  claim  that  this  condition  Is  a  sufficient  condition 
for  stability  when  taken  together  with  v  >  0  and  eq.  (62). 
For  suppose  the  first  product  is  positive  at  the  value 
of  z  In  question.   Then,  upon  expanding  the  square  In  the 
last  term  of  (67),  we  obtain,  after  collecting  like  terms, 
a  coefficient  of  the  first  product  which  is  given  by 

CO  00   I    I 

(69) 


^-1 


o,v'2 


Po  "°'^  2 


Ia-P-|  /  U-^ 


dfi. 
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o  /A 


v=z 


which  is  less  than 


/   (-X,v)e^^/   (^fo,v)x<^- 


2 


+  + 


|a  P"l  /   M- 


Po  A 


dM- 


V=2 


if  eq.  (68)  is  satisfied.   Therefore,  under  the  hypothesis 
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of  positlvity  of  the  first  product  and  the  validity  of  eq.  (68), 
we  have,  at  v  =  z. 


(70)        Re(D(2;A^))    <  - z^A^   ^ 

^o   I  J 


o,v1 


d[x 


o     /2 


(^^^^Cv^S   ^^^ 


+  + 
-|a  p 
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\'i^-^    aixY 


o 


'A       J 


p       '      "^  o 
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/  M 


'f 


'o    /a 


dkL 


^T-'-^'^''4^iv'oM 


^o    /Z 


-   A^  -^    (Re   T)         (H^f^      )      d[i(        <    0 


^ 


o,v'2 


00 

If  /   (m-  f  „/p„).  dM-  =  0,  but  our  Inequality,  given  by 

O , V/   O  A 

eq.  (68),  Is  satisfied,  the  remaining  two  terms  have  the  same 

2 
sign.   From  Im  D(z;A  )  =  0  we  obtain 


(71) 


B^ 
0  =  A^  ^  (Re  T(z)) 


^  (fo,v)s  dt^ 


^4-^^'=44v''-t 
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Po  '°'^/E 
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together  with 
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(72)  Re   D(z;a2)    J,^^   ^^    ^^  _y^  _j,^  ^^^   ""   P^  ""   P^  V^'^°^^^^''i 


-   A     — ^  z      a   p       si    ' — 


2,2    -^o      2 
,    —      ,  ,     ,       -  TT    A       Z 

°  1  ^^vpo  ;a\         Po 


00 


d[^ 


a 


Po      °'VZ 
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X 


^    (fo,v)E    ^^ 


However,    since   /     |J.    (f        )      dM-     and      /    (-L^  f        )      diJ.     have 

o      °'^  ^  o   Po    ' 

the  same  sign,  from  eq.  (71)  we  Infer  that  Re  T(z)  and 

o    2 

/     Q       Q       Q       O         -^  X  /-\  ) 

z  -V  -A  (c   +  ^  +  ^  i  2^^o  v^2M   '^^^^  opposing  sign 

^  Po     Pq   l-L      '  J 

2 
and,  therefore,  from  eq.  (72)  we  conclude  that  Re  D(z;A  )  <  0. 

2  2 

If  A  =0,  then  on  the  real  line  D(z;A  )  becomes 


(73) 


D(z;0)  =  (z2-v2)[l+  lal^-  z^I  (— ^-^^ 


-  TTIZ 


a  f 


o,v 


dii.; 


o  /Z 


o 


so  that  the  only  real  values  of  z  different  than  zero  for 
which  D(z;0)  Is  real  are  given  by  roots  of 


(74) 


a  f 


o ,  V 
'o    IZ 


d|-L 


0  . 


V  =  Z 
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However,  if  eq.  (7^j  Is  satisfied  the  Inequality  given  by 

eq.  (68)  Is  not. 

2 
Therefore,  v  >  0  together  with  the  Inequalities  given 

by  eq.  (62)  and  (68)  are  sufficient  to  guarantee  exponential 

stability. 

We  note  that  In  the  test  prescribed  In  eq.  (68)  the 
kinetic  properties  of  the  equilibrium  configuration  are  pre- 
served, and  as  a  result  this  sufficiency  condition  Is  very 
sensitive  to  the  kinetic  velocity  profile  of  the  distribution 
functions  In  \i-v   space. 

We  now  describe  certain  classes  of  functions  which 
satisfy  this  criterion. 

+       O  p 

If  f~(M.,v  )  Is  monotone  down  In  v   for  every  |-l,  then  we 
claim  that  the  Inequality  given  by  eq.  (68)  Is  always  satis- 
fied.  This  follows  from  an  expansion  of  both  sides  of  the 
Inequality.   On  the  left  we  have 
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+     2 
from  the  monotonlclty  of  f^(l^,v  ).   On  the  right  we  obtain 
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(76) 
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using  the  monotonlcity  again.   Using  Schwartz's  inequality 

I  ^    ^il/2       I  t    +,1/2 
with  the  functions   f    P    '   and  I-l  f    P^  '  ,  we  have 

that  the  right  hand  side  of  the  inequality  given  by  eq.  (76) 


is  greater  than  that  given  by  eq.  (75)-   This  completes  the 

proof. 

-     2  2 

Suppose  now  that  f  {[i,v    )  are  monotone  down  in  v   for 

+  + 

|v|  <  v~(m.)  but  are  identically  zero  for  |v|  >  v  (ti.)  (i.e., 

v-(m.)  Is  a  cutoff  velocity  depending  upon  I-l).   Then  if  there 

+ 
exists  a  maximal  v   for  v  {[.i)   which  is  finite,  which  we 

denote  as  v~,  then  we  observe  that  at  z  ^  v  ,  the  inequality 

given  by  eq.  (68)  fails  and  if  Re(D(v  ,A  ))  >  0,  a  contour 

diagram,  shown  in  Fig.  9  (or  something  qualitatively 

equivalent  to  it)  might  occur.   We  make  the  claim  that  this 

cannot  happen.   That  is,  functions  which  are  monotone  down 

2 

in  V  ,  modulo  a  cutoff  velocity,  are  also  exponentially  stable 

The  proof  is  simple.  Suppose  these  functions  were  not 
stable  exponentially.  Then  in  a  small  half  neighborhood  of 
V   given  by  v  <  v  ,  we  must  have  Im  D(v,A  )  >  0.   That  is 
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to  say,  the  contour  diagram  must  produce  a  crossing  of  the 
real  positive  axis  and  not  just  touch  it.   However,  this 
implies  that  in  a  neighborhood  of  this  line  in  the  complex 
upper  half  plane,  the  imaginary  part  of  D(z;7\  )  must  also  be 

positive,  by  the  continuity  of  D(z;A^).   Consider,  however, 

+     2 
the  sequences   g~(M-,v  )   defined  by 


''n 


(77) 


gn(^^'V^) 
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o,v^    n 
±  +    +2 

V  f~(l-L,V    ) 

n  o^  '  n  ' 


;v^-v^) 
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f:(^^.v  ) 


■  o 


V. 


if  |v|  <  v-(n) 


+       +  +       ± 

where  v  (^l)  <  v"(^i)  and   11m  v^(^)  =  v^(^i).   We  have  that 


(78) 
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n- 


•'n 
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lim  g 
n-^  00 


n,v 


+ 
"  o,v 


t  2 

The  ff  are  monotone  down  in  v   for  all  values  of  \i   and, 
^n 

therefore,  satisfy  our  sufficient  condition  for  stability. 
However,  for  every  z  belonging  to  the  upper  half  plane,  the 


mapping  by  D(z:A  )  depends  continuously  on  the  functions 

+      4- 

f~   .   Since  the  image  of  the  boundary  of  the  half  plane,  with 

+  +  2 

g~  replacing  f~,  under  D(z;A  )  does  not  cross  the  positive 


r 


See  appendix. 
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real  axis,  we  have,  for  fixed  z  and  sufficiently  large  n 
that  near  v~,  the  Image  of  z  Is  In  the  lower  half  plane. 
Thus,  In  the  limit,  the  image  must  be  in  the  lower  half  plane. 
This  completes  the  proof. 

We  note,  before  leaving  this  subject,  that  the  method 

of  proof  also  yields  the  exponential  stability  of  those 

2 
functions  which  are  only  monotonically  non-increasing  in  v 

for  every  value  of  [x,  provided  they  satisfy  the  other  hypo- 
thesis we  have  imposed. 

We  conclude  this  section  with  a  treatment  of  the  problem 
excluded  in  the  previous  discussion  of  stability.   That  is, 
we  assume  the  symmetry  of  the  problem  to  be  such  that  all 
spaclal  dependence  of  the  solution  is  confined  to  the 
variables  in  the  plane  orthogonal  to  B  • 

In  this  situation,  our  equations  become 

+ 


S    +   ^ 


(79)       ^(f--B^B,)=0 


(80)      ^  -^  (sT  -^2  +  ST  --j)  =  ° 

2         5 


(81)       ^^2,^-^ 


^^2,3  ,   2  ^^1       ,  ^l         S 
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M-dndv(f^  +  f  )  =  0, 


(82)       -^^  +   c"  ^— t-  +  ^  ^ 

Equations  (79)  and  (8l)  are  easily  integrated  and 
yield 
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B^^^{t,x)    =  62^3(0, x) 


(83)  + 

+  f  + 

f  {t,x;[i,v)    -  g^B^(t,x)  +  P  (0,x;^i,v) 

o 

+ 
where  F~  is  a  linear  combination  of  the  Initial  conditions 

and  Is  time -Independent . 

Substituting  the  second  of  eq.  (83)  into  eq.  (82)  we 

obtain 


(8*)      tI^  Mc^  +  p-)  ^  =  -  33^  G(0,x) 


where 


b2  r 

G(0,x)  =-^ 


M.d[i.dv[P"''(i-L,v)  +  P  ([i.,v)] 


Assuming  that  our  functions  have  two  continuous  derivatives, 
we  may  eliminate  rp  and  r^  to  obtain  a  wave  equation  for 
B.  with  wave  speed  \ —  (B  /M-  +  2p°( 


^^B.      p    p°    ;x2     .2        .2     .2       . 
(85)      -^  -  (c^  +  ^)(-^  +  ^)B.  =  (-\  +  -\)  G(0,x), 

ht  ^o      bx^        bx^  Sx^    hx^ 

25         23 

The  solution  to  this  equation  and  its  properties  are  well 

17 
known.     It  suffices  for  us  to  note  that  the  wave  speed  is 

real  and  that  no  additional  criteria  for  exponential 

stability  are  introduced  in  this  singular  case. 

The  appearance  of  this  wave  speed  will  be  accounted  for 

in  the  next  section. 
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V.   Some  Asymptotic  Results  Pertaining  to  Wave  Propagation 

We  recall  that  In  the  treatment  of  propagation  orthog- 
onal to  B  ,  It  was  observed  that  signals  propagated  with 

f  1  °  2        o  7 1/2 
speed  \-^   {B  /yi     +   2pj_)>    .   Since  a  cursory  Inspection  of 

o^^o 
D(z;A  )  does  not  predict  such  a  wave  speed,  we  now  will 

show  that  the  result  Is  obtained  from  this  function  by  an 

2 
asymptotic  expansion  of  A  . 

Defining  the  wave  speed  o   =   z   cos  0  =  '^^/ 1  k  |  ,  we  note 

that  non-zero  a  and  0   — >  7r/2  Implies  that  |z|  — >  co.   Thus, 

we  allow  A   to  become  large  (0  — >  Tr/2 )  and  look  for  non- 

zero  wave  speeds  which  yield  zeros  of  D(z;A  ).   Since  we 

have  just  seen  that  the  existence  of  such  a  propagation 

speed  Implies  |z|  — >  oo,  we  may  employ  the  asymptotic  expan- 

1"        2 
slon  of  D(z;A  )  for  large  z  and  obtain. 


D(z;A^) 


+ 


J 


n; 


o 


o/A 


o 

2  f    2   P^x  ,   2, 
V  +(c  +  — )  tan  ( 


l-Ldi-Ldv  f. 


Factoring  l/cos  9   and  retaining  only  the  dominant  part  of 
the  expression  we  see  that 


D(z;A^)' 


1 


2a 
cos  9 


_o 

a   -  (c  +  -— )  sm  i 


t 


See  appendix. 
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and  this  is  the  desired  resijlt. 

Some  further  asymptotic  results  may  be  obtained  If  It 
Is  assumed  that  the  ratio  of  Ion  and  electron  thermal  speeds 
Is  much  smaller  than  unity  (i.e.,  v,,/v,,  «  l)  and  we 
Investigate  the  existence  of  zeros  of  D(z;A  )  which  are 


+ 


large  with  respect  to  v,,  but  small  when  compared  with  v 


th 


+ 


We  will  assume  that  f  (|J-,v)  are  of  the  form 


th' 


f-(^.v)  =  -^ 


+ 


1  2 

2  ^ 

('^th)       (^th) 


B  l^ 
o 


+ 


and  define  the  quantities  ^   as 


+ 

r  = 


z 


thereby  being  led  to  soli; 
Defining  the  dlmens'- 


:  with  II  I  »  1  and  ||"|  «  1 
J 3  variables 


y 


B^H 


V 


th 


X  = 


V 


■th 


we  obtain 


r 


+ 
f - 

k   o.v  -  -, 
v-z 


o  /  2  xk-1 


J 


B 


k  (^th) 


r  k  -,   1  2. 

y  g^iy.  2X  ) 


J 


xA 


T" 


dydx 


Therefore,  in  terms  of  these  dlmensionless  variables,  we  find 


-  ^h  - 


(86) 


D(z;^^)  =   (l+|a|^)  -  z^  (i|- 


V 


th 


5^  dydx)^ 


^^-(v^-^^(c^-^)-^^(PovL[S^^^-),^ 


-  2  p^  (|a|p^)  <)  ( 


O' 


y  3^  dydx)^ 


•  +i 


We  now  make  use  of  the  assumptions  | ^  |  »  1  and 
U  I  «  1  to  obtain  the  following  expansions  of  the  Integrals 
In  eq.  (86) : 


r  k  + 
.+ 


(87) 


dxdy 


r 


x-C  ^ 


r  k  - 


J 


x-^- 


dxdy 


1_ 
2 

+   (  J 

k 


y  g  dxdy  + 


'+ 


2  k   ,  ,    ^/  1  N 
X  y  g  dxdy  +  O(^rj-) 


7T> 
+  J 


J 


y  g  2  dxdy  +  ttII 

X 


y  gxx(y.o)dy 


+  0(^^) 


where  we  have  assumed,  for  convenience,  two  derivatives  for 
g~- 

Substituting  the  expansions  of  eq.  (87)  Into  eq.  (86) 
we  obtain  the  asymptotic  result 


.2.. 2,  2 


(88) 


D(z;a2)  ^+1(1+   \a\     )\P 


th' 


f      v-^^-x-^cc-^  ^  p;/p„) 


I 


+ 


V 


th 


'V. 


-  \2 


P,-  ^o|_  + 


o    W 


th 


2  -   ^  ^  (    e2  A""  . 

y  g  2  ^^^y    -  ^+  p-  ( 


+ 


X 


o 


a|po)~(|yg  2^^^y)  ^i 

X 
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In  eq.  (88),  only  the  dominant  terms  have  been  retained.   It 

2  t 
has  been  assumed  that  (P^v^^)   are  of  the  same  order  of 

magnitude,  and.  In  addition.  Integrals  of  Ion  and  electron 

functions  which  appear  are  of  the  order  of  unity.   The  first 

of  these  assumptions,  together  with  charge  neutrality  and 

the  hypothesis  v^,/v~  «  1  allows  us  to  assume  that 


m 


/m"*"  «  1  so  that  we  may  take  |a|   =0  and  |a|   =  1, 


together  with  p   =  p   for  our  purposes. 

Writing  ^  =   i    cos  6,    we  find  that  to  lowest  order  the 
zeros  of  D(z;A  )  are  given  by 

T^Kv   cos  0+  (c  +  - — )  sm 


^2  ^]4J^ ^ 


—  ^  (/  yg  2  ^y^^) 


Po 


+  ; X 

1  1  -  V  -^  (  /  ys~2  dydx) 
L       Pq      ^ 

2   —   + 
If  A  Pr/Po  '^'^  "*-  (^-^-^  ^  ^s  ^°^  near  t/2)  ,    the  denominators 

may  be  approximated  by  unity,  so  that  in  order  that  |^  |  »  1 
we  must  have 

1  (f^2/,,      ,   Ox  ,   o   .  2„     o    2, 
—  [(ByM-Q  +  P^)  +  Pi  sm  9  -  p^^  cos  ( 


Po 


+2 
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»  1 


2 

Assuming  that  p"^  v^u  is  of  the  order  of  the  fluid  pressures, 
^  "^  o  th 

another  way  of  expressing  the  result  Is  that  for  the  waves 
to  exist  It  Is  necessary  that  the  ratio  of  the  magnetic  to 
the  fluid  pressures  be  much  larger  than  unity.   This,  of 
course,  corresponds  to  small  p  In  Pigs.  7  and  8. 


VI.   Fundamental  Solution  for  the  One -Dimensional  Problem 

We  proceed  now  to  find  a  fundamental  solution  for  the 
one-dlmenslonal  problem  associated  with  equations  (11) -(13)' 
It  Is  assumed  that  the  Independent  space  variable,  s,  lies 
on  an  axis  In  the  x-y  plane  making  an  angle  9   with  the 
magnetic  field  line.   Under  these  hypotheses,  the  equations 
become 


B  ,  .  a    +    ^       +        + 


(89)  0  =    (^+  V   cos   0  ^)f-  +    ([nJa"-   \^)f~^^ 


V 


+  t 

-   V   ^  [cos    0  ^  B,   +  /  sm   e   I-  r^   +  a-f^^^   cos    9  ^     v' 
o   j  o 


(90) 


0   =  ^  B^  +    sm   9   ^  r^ 


°   =^^^2,5    -   ^°"    ^   ^''2,3 


-   ^7 


(90) 

0  =  ^  -  V  cos  0  ^  B^  . 


We  define,  as  before,  the  vector 


(91)      u(t)  =  (f^,  f  ,  B^,  B^,  v^,    B^,  r^) 


with  Initial  data 


_\. 


+ 


(92)      u(0)  =  (g^,  g",  b^,  b^,  Rg.  b^,  R^)6(s) 

Since  the  problem  for  9   =  it/2,   has  been  solved  explicitly 
In  a  preceding  section,  we  assume   9  /   7r/2  and  divide  our 
equations  by  cos  0.   Introducing  the  new  variable  t  cos  0  =  t, 
we  Fourier  transform  the  equations  In  s  and  take  a  Laplace 
transform  In  t  to  obtain 


(93)       (p+  lkv)f"  +^^([T]]a"-  ^i)fo^^]  ikB^ 

+ 

+  IkA  3^  ?2  +  "'^o,v  ik<v^)  =  S"(^^,v) 


(94) 


pB-,  +  Alkr   =  b-. 


^^2,3  -  ^^"^2,3  ^  ^2,3 
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B 


(95) 


pr^  -  ikv  Bg  +  ike  'KB-^  +   IkA  ^-^Xl-i-/  ^   Ro 


pr^  -  Ikv  B^  =  R^ 


where,  as  before,  p  is  the  Laplace  transform  variable,  k 
is  the  Fourier  transform  variable,  and  A  =  tan  9. 

We  observe  that  eq.  (93) -(95)  are,  after  defining 
p/ik  =  -z,  exactly  eq.  (^3)-(5^)  with  a  slightly  modified 
right  hand  side.   We  readily  obtain,  therefore,  eq.  (46) 
with  a  modified  right  hand  side  given  by 


(96) 


+  .  - 


^  (/  n  ^J:^  dii-dv,  /  V 
ik  ^  -^    v-z      '  •' 


2  gVPo-s"/Po 
v-z 


Writing  the  matrix     A     as 


dn-dv,    b-,  , 


^2  j>    ^2'    ^3'    ^3) 


z' 


(97) 


M 


0 


V 


0 


-z      -1 

2 

-V       -z 

J 


as   in  eq.    (48'),   we  have 

r 

]y~i 
-1 


(98) 


0 


0 


2      2         2      2 

Z     -V  Z     -V 

2 
V  z 

~2      2  2      2 

z     -V  z    -V 


y 
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.-1  _  _  ^ 

zD 


(99)  M'^   =   -  ^    >^ 


'\ 


-z(z2-4)T      -z{z^-4)l^{\a\f^^^)^      -(z2-v2)(ilj  v^[i]-za[i 


z  !^  Itt]    -.r.2_„2_^2   ^  (^^]^         _^^2_^2^^^^|         v^A^Ill-zAUl 


Po 


2^2  g2 

T        -zA^  -^  I    f  |a|f^   ^J.         (z^-v^)T  -v"AT      -zAT 


^   Bq  ^2   ^o    _    ,,     ,  .         X  /_2    ..2^^  .2 


Pq  Po      ^^''     '     °'V'^ 

2  2 

zA  ^  T  zA^l(la|f         )^  ^"(^    ""^    ^^  D+v^T      -zT 

p  p  IJ.^  '       '      OjV^A  A 

2  2 


o  "^  o 


V 


where 


o  4-  f 

(  B^  M-  o 


"I'itll^-^V'Tf'A 


2  2,2    2 

U)■^  =        V  +        A        C 


We  shall  first  treat  propagation  parallel  to  the  magnetic 

—  2   2 

field  vector,  B  .   In  this  situation  A  =  0  and  D(z;0)  =  (z  -v  )T. 

As  a  result,  M   is  given  in  the  relatively  simple  form 
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x 


(100) 


M 


_-l 


+ 


(l-IPo)"t^JTf^J  ^,j 


^ 


0 


0 


0 


0 


0 


0 


1 

T 
0 

0 

0 


llli 


zl 

\j 

1 
z 

0 

0 

0 

-z 

- 

2      2 

Z     -V 

?      P 

Z     -V 

0 

p 

-z 

^      2 

2      2 

Z  -V      Z  -V 


V 


J 


If  Q  is  one  of  our  transformed  dependent  variables, 
then  q(T,s)  is  given  by 


(101)  q(T,s)    =  Re  -k^   J    kdk 

2^'  J 


+oo4-ia 


ikfs-zx)    ^  f      ,  \  , 
e      ^  ^    Q(z,k)dz 


-a>+xa 

where  a  is  chosen  sufficiently  large  so  that  all  residues  of 
q(z;kj  have  imaginary  part  less  than  a.   We  shall  assume 
that  T  has  no  zeros  in  the  upper  half  plane  Im  z  >  0. 

We  first  observe  that  for  all  finite  values  of  A, 
we  have 


oo-fiCT 
r 


B^(t,s) 


=  Rt 


2.2i 


(102) 


dk 

o        -o>fia 


^ik(s-ZT)    _|Z_^    (^^_^        ) 
Z     -V  ^  '^ 


I 


r^('^>s) 


R£ 


2Tr~i 


dk 


_oo+ia 

ik(s-ZT)      dz 


2-^   (V   b    -zR    ) 

Z    -V 


O  -oo+xc 
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Using  the  definition  of  the  5-function 


18 


+  00 


(103) 


5(x)  =  Re 


IT 


ikx  -,,     1 
e     dk  =  2? 


ikx  -, 
e    dK 


we  Immediately  obtain 


B^(t,s)  =  2"  t)^[5(s+VT)  +  5(s-vt)] 


l!3 

2  V 


[5(s+vt)  -  &(s-vt)] 


(10^) 


r^(T,s)  =  p-R^[5(s+VT)  +  5(s-vt)] 
-  p-  vb^[5(s+VT)  -  &(s-vt)] 


which  are  simply  f-undamental  solutions  to  the  wave  equation 
with  vjave  speed  v.   Similarly,  we  obtain  identical  results 
for  BoC^jS)  and  r^Cx^s)  by  simply  changing  the  subscripts  in 
eq.  (104-)  from  3  to  2  for  this  special  angle  of  propagation. 
For  B-,(t,s)  we  easily  obtain  from  eq.  (90 ) 

(105)     B^(t,s)  =  b^5(s)  . 

The  computations  for-^^U-^  and  <^v  \  are  m.ore  subtle. 
We  have 


52  - 


(106)  (\iy  -  Re   -^ 

2ir  1 


00  00+ia 

r 


dk 


dz   e 


Ik(s-ZT) 


\x  ^  "'"^     di^dv 
v-z 


o        -oo+ia 


+ 


+ 


(   a  p    )  f 


T 


^i^ ^  +   ^  j  I 

v-z  z     I 


Here,  we  are  unable,  in  general,  to  continue  the  integrand 
analytically  into  the  lower  z  half  plane.   We  note,  however, 

that  the  k  integration  may  be  continued  bo  the  negative 

19 
k-axis  without  contributing  to  the  integral,    and,  there- 
fore, the  range  of  the  integration  may  be  taken  over  the 
entire  real  k-axis.   Utilizing  the  fact  that  T  has  no  zeros 
in  Im  z  >  0,  we  displace  the  contour  of  integration  in  the 
z-plane  to  the  real  line  and  obtain,  changing  the  order  of 
integration 

+ 

f  r  (   \ry\n       ) 

(107)    <n(T,s)>=i 

f 


T 


M'gy(i^,   Y^^^  + 


(  |ct  I  P  ) 


^Z 


TT 


Im 


I  (  Q^V)  iv^(^)  ' 
V^  P^  ^A^^  ^p^^A 


+  b. 


s 
z=  — 

T 


O 
-5-+  I  2 


^fo,v)z 


B 


\J- 


V 


+  P^  (1)  i^i 

IT  ^S  ^      •- 


6(s) 


Similarly, 


(108)    <(v^(t,s)) 


if  (I^Ipq)' 


T 


+ 


IT 


Im 


I(v''(g/pj 


o^A 


T 


Z  =  S/  T 


/' 


-b- 


[n]5(.)  +  i  (f)lm 


riM.'^o,>o'A 


z=s/t  -^ 
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We  note  that  in  eq.  (107)  if  we  assume  that  f  (|J-,v)  is  even 
in  V  with  a  partial  derivative  in  v  which  goes  to  zero  as 
vp  (p  >  1),  then  the  principle  value  notation  is  not 

necessary  since  Im^l]  consists  of  a  sum  of  terms,  all  of 

~  k    t 
which  contain  a  factor  of  the  form  /  [^   di-i  f    (l-i-,v) 


(k  =  0,1,2) 


o       °'^ 


v=s/t 


The  behavior  of  B  and  r  for  the  various  ranges  of  s 

and  T  are  easily  found  from  the  6-functions.   Specifically, 

Bo  V  and  r^  -,  are  zero  except  on  the  wave  fronts  s/t  =  +  v  , 
2,5       2,p  -  ^  -       ' 

along  which  an  undamped  disturbance  is  seen  to  propagate. 

B.  is,  of  course,  constant  in  time. 

Inspection  of  eq.  (IO7)  and  (I08)  shows  that,  in  general, 

on  every  wave  front  s/t  =  constant,  there  is  a  signal  which 

is  observed  to  damp  in  time  like  t   .   Fixing  the  time  and 

considering  the  behavior  of  the  signal  as  experienced  at 

points  in  physical  space  at  large  distances  from  the  initial 

disturbance,  '\1^)'  and  <^v  /   are  seen  to  decay  in  distance  like 

/  1-L  (|    )(|j.,  s/T)d!i.  r      ;  speaking,  so  that  the  effect 

o     o,v 

of  this  disturbance  is  to  be  attritubed  to  the  fast  particles, 

as  was  to  be  expected.   On  the  other  hand,  if  we  consider  the 

effect  of  the  signal  at  any  single  position  after  a  large 

time,  we  observe  a  much  slower  dam.ping,  with  asymptotic 

behavior  in  time  dominated  by  the  slow  particles.   In  fact 


(109)      <^L>  ~  ^    ug^(n,0)dM- 
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(109) 


<^'>'^(l+|aL)  J2    N 


Z'  T 


^\ 


Vt(^-^/^) 


di-L 


o 


A 


+ 


alp  ) 


&l^\    a,) 


O 


A 


Assuming  B.  and  <\V  ^  are  known  functions  of  z,  we 
easily  solve  for  f   In  eq.  (93)  to  obtain 


(110)      f-  =-^  f^  -  F-(^x,v)B3_  -a  f-^^<v  > 


with 


+ 


+ 


+ 


o 


P"(^^,v)  =   ([Ti]a--  \^)^l^^    -   V  g- 


Uslng  the  previous  methods,  and  accounting  for  the 
singularity  on  the  real  line  at  z  =  v,  we  find 


(111) 


+ 


■+       +  ^       +  + 
g~  -  Re(F~B  (v)  +  a  f 


o,v 


v^(v)f)^6(s-vt) 


+ 


\   P-V.  ^  Im(pV(f)  +  :-i\^^   <v"(|)>  ) 


where  q*(z)  =  lkq(z,k). 

In  the  same  way  we  obtain  the  representation  for 

A  >  0.   The  major  distinction  between  the  two  cases  is  that 

2 
for  A  >  0,  the  effect  of  a  disturbance  in  Bp  and  v^   is  felt 

+ 

by  f  and  B-,  and  vice-versa.   That  is  to  say,  Bp  and  v^ 

remain  coupled  to  these  variables  when  the  wave  vector  is 
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not  parallel  to  B  •   Also,  If  D(z;A  )  does  not  have  real 
zeros  in  the  z  variable,  undamped  wave  propagation  Is  only- 
observed  in  the  solutions  for  r^  and  B^,  which  are,  of 
course,  still  independent  of  disturbances  in  the  other 
quantities . 

Rather  than  writing  down  the  solutions,  which  are  very 
lengthy,  we  present  the  recipe  which  is  valid  for  all  of 
the  quantities  ^M-),  <(v  )> ,  B-,  ,  Bp ,  and  r^.      One  need  only 
Identify  the  appropriate  elements  of  the  matrix  in  eq.  (99)' 

Let  q(t,s)  be  one  of  the  five  variables  mentioned  above 

p 

If  D(z;7^    )    has  no   zeros   for  real   z,    then   q(z;k)    is    of  the 
form 


(112)  a(z,k)    =  jk|A^(z,A) 

+   A^(z,A) 


/• 


M-  -^^—  dl-Ldv 
v-z 


^v^(^p    )                       A.(z;A)                A,(z,A)  ^ 

— TT^  d^dv   +   ^ b^  +  ^4 b^   +   A3(z,A)rJ 


where  A  (z,A)  is  analytic  in  z  for  Im  z  >  0  and  bounded  for 
Im  z  >  0.   When  q(t;s)  =  rp(t,s),  A^  and  Aj,  contain  a 
factor  z,  so  that  the  exhibited  singularity  at  the  origin  does 
not,  in  this  special  case,  exist. 
Thus,  we  obtain 

(113)    q(T,s)  =:^   Im|A^(z,A)l^(g^)  +  A^ (z , A)lv^ (g/p^)^ 

A^(z,A)b^+A^(z,A)b2 

z=s/t 

+  A^(0,A)b2|6(s) 
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+  P.V. 


2 I  +  A  (z,A)R2V      +  Re<A^(0,A)b^ 


+ 
The  representations  for  f  are  obtained  exactly  as 

before,  with  the  trivial  exception  of  the  addition  of  a 

term  due  to  the  fact  that  r^Ct^s)  does  not  separate 

with  Bp(t,s)  from  the  system. 
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Appendix 

¥e  will  be  concerned,  in  this  section,  with  properties 
of  functlonals  of  the  form 


+  00         oo 


(A.l) 


-  k 

l^   J 


2n, 
V   ( 


00        O 


^^r       dn)  ^^ 

o,v    '  v-z 


when  z  is  a  complex  number  in  the  half  plane  Im  z  >  0, 
n  =  0,1  and  k  =  0,1,2. 

If  n  =  1 ,  we  may  reduce  the  above  integral  to  the  case 
n  =  0  by  adding  and  subtracting  z  twice.   Thus, 


+  00         00 


+  00 


(A. 2) 


J      J 

-00      O 


p.kf      an)     ^^ 

o,v    ''  v-z 


v(v-z)+vz 


v-z 


J 

—  00 


\i^t         diJ.)dv 
o,v    ^ 


+  00         00 

r        ^ 

vdv 
1 

-00      O 


M-^f    d^L  +  z( 
o,v        ^ 


+  00       00 

dv 

-00     O 


^^fo,v  ^^)  +  ^^ 


+  00         00 

dv    k„     , 
—=-^       1^  f    dv 
v-z      o,v 

'-00      O 


After  an  integration  by  parts,  we  obtain 


(A. 5) 


+  00  ^ 

2   ^"^ 
V     v-z 


+  00  00 

r  r 


J 

—  00 


J 


k  2 

^i  f      d^i  =  z 

o,v 


o 


J 


dv 
v-z 
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When  k  =  0  or  1 ,  the  last  term  is  a  mass  density  or  a 
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perpendicular  pressure  divided  by  the  magnetic  field. 

Expanding  the  denominator  of  (A.l)  In  Inverse  powers 
of  7.,   we  have 


(A. 4) 
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so  that  for  large  values  of  z,  we  obtain 
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(A. 5) 
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In  particular.  If  f^(ii.,v)  Is  even  In  v,  the  first  term  In 
the  expansion  vanishes,  so  that 


(A. 6) 
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¥e  now  show  that  If  ]/  lo-  f^"{  dM-j  is  a  convergent  sequence 


o,v 


In  L   (p  >  1)  then  the  sequence  of  functlonals  lY   obtained 

P  , ,  K 

(m)  ^ 

by  replacing  f^  ^  In  (A.l)  converges.   In  view  of  (A. 3)  It 
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suffices  to  prove  this  for  n  =  0.   Let  z  =  x+ly.   Then 


(A. 7) 
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Since  I  /  M-  di-L  f^"^-'   is  Cauchy  In  L  ,    this  completes  the  proof, 
If  we  assume,  once  more,  that  f  (ix,v)  Is  even  in  v,  then 
Re  I  ,  is  even  and  Im  I  ,  is  odd  in  Re  z.   Again,  referring 
to  (A. 3)^  we  need  only  prove  this  for  n  =  0. 


(A. 8) 
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Suppose  in  (A .8)  x  is  replaced  by  -x.   Changing  variables 
(v  =  -v  ) ,  we  have  |v-z|  =  |v  +z| ,  (v-x)  simply  changes  sign, 
so  that,  replacing  the  dummy  variable  v  by  v" ,  the  real  part 
of  the  integral  is  unchanged  but  the  Imaginary  part  is 
negated . 


CO 
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We  now  show  that  f^  =  f^d-^BQ  +  ^v    ).      Then  Pj_  =  p^^  . 
Prom  the  definitions,  we  have 
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Letting  |J.B  +  ^  v   =  |  and  changing  the  orders  of  integration. 
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we  obtain 
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Q.E.D. 
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Domain  £) 


z-plane 


Figure  I 


Image  of  £) 


C-plane 


ImC  =  0 


Figure  2 
Illustrating  a  One-One  Mapping  of  4)  by  f  (z) 
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Values  of  C  which  ore 
assumed  once  by  f(z) 
m  JD 


C-plane 


Values  of  ^  which  are 
assumed  twice  by 
f(z)  in  £)  


••lnn^=0 


Figure  3 
Illustrating  a  Mapping  of  £)  by  f(2)  Which  is   not  One-One 


Values  of  C  which  ore 
assumed  once  by  f(z) 
in   £)  


Values  of  C  which 
are  assumed  twice 
by    f(z)  in   £) 


^-  plone 


*-Im{=o 


Figure  4 
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Domain  of  pressures  for  fixed 
magnetic  field  strength  which 
permit  stability  (0<r)<G) 


Pn--T   p!-    (^  +  2P, 
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Figure  5 


BV/Xc 


Domain  of  pressures  for  fixed  magnetic  field 
strength  which  permit  stability  (7;>6) 
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(^-2) 


Figure  6 
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Range  of  P||/P^  which  permit 
stability  for  given  values  of  /3. 
(0<M/6<i)    \  \ 


Figure  7 


97/6 


P  /P 


R/Px='  +  ^ 


Range  of  P(|/Pl  which  permit  stability  for  given 
values  of  /S.   (7^  >6) 


2.(/rr4^-3) 


Figure  8 

Illustrating  Existence  of  Maximal  /3  Beyond   Which  Stability  is 
Impossible  if    7^  >6 
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,Re£  =  0 


Values  assumed 
twice  by  D{z;X^) 
in   JO 


•"Tm^  =  0 


Values  assumed  once  by 
D{z;  X^)  in    JO 


Figure  9 

Illustrating  an  Impossible  Mapping  Induced  by  Cut-Off  Functions 

Monotone  in   v2 
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